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Abstract 

In the present work, ternary Special Quasirandom Structures (SQSs) for a fee solid solution phase 
are generated at different compositions, xa = xb = %c = § and xa = |, xb = xc = \, whose 
correlation functions are satisfactorily close to those of a random fee solution. The generated 
SQSs are used to calculate the mixing enthalpy of the fee phase in the Ca-Sr-Yb system. It 
is observed that first-principles calculations of all the binary and ternary SQSs in the Ca-Sr-Yb 
system exhibit very small local relaxation. It is concluded that the fee ternary SQSs can provide 
valuable information about the mixing behavior of the fee ternary solid solution phase. The SQSs 
presented in this work can be widely used to study the behavior of ternary fee solid solutions. 
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I. INTRODUCTION 



In CALPHAD (CALculation of PHAse Diagrams) modeling 1,2 , thermodynamic properties 
of a solution phase in a ternary — or higher-order — system are usually obtained through 
combining those of its constitutive binary systems and occasionally with ternary interaction 
parameters. Since the most dominant interatomic reaction in a multicomponent system is 
that of the binaries, accurate thermodynamic descriptions which are capable of reproducing 
the characteristics of binary solution phases are prerequisites to a successful multicomponent 
thermodynamic modeling. In this regard, considerable efforts have been made to develop 
models for combining thermodynamic descriptions for binary solution phases to be used in 
higher-order systems 3 ' 4 ' 5 ' 6 ' 7,8 . 

In typical thermodynamic modelings, ternary interaction parameters for solid solution 
phases are not introduced since their thermochemical data, such as mixing enthalpy, which 
are needed in the parameter evaluation are very difficult to measure. In fact, obtaining 
accurate thermochemical data for solid solution phases is challenging even for binaries due 
to sluggish kinetics at low temperatures so that it is hard to reach a complete thermodynamic 
equilibrium. Furthermore, the existence of intermediate phases narrows the experimentally 
accessible composition range for thermochemical properties of solid solution phases. As 
the number of elements increases in a multicomponent system, the complexity of acquiring 
reliable data also increases. Consequently, interaction parameters for the excess Gibbs energy 
of binary solid solution phases are usually evaluated only from phase equilibrium data and 
those of ternary are usually omitted due to the lack of data. 

Fortunately, the dearth of experimental data for solid solution phases can be overcome by 
atomic level calculations, such as first-principles calculations. Specially designed supercells, 
so-called Special Quasirandom Structures (SQS) suggested by Zunger et al. 9 , whose lattice 
sites are occupied by constituent atoms so as to reproduce the average correlation functions 
of a completely random solution, are one of the representative methods. SQSs can be 
completely relaxed within density functional theory (DFT) codes to consider the effect of 
local relaxation and can also be applied to any other system by changing the atoms because 
they are structural templates. A limitation is that SQSs of a unit cell size that is manageable 
with present DFT codes can only be obtained at certain compositions, e.g. £=0.25, 0.5, 
and 0.75 in the substitutional A\_ X B X binary alloys. Nevertheless, first-principles study of 
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SQSs at those three compositions permit accurate predictions of various properties of solid 
solutions. It has already been successfully applied to calculate thermodynamic properties of 
binary solid solution phases for fee, bec and hep phases 10 ' 11,12 . 

Two ternary fee SQSs in an A-B-C system, which can be applied to four different com- 
positions at x A = x B = x c = |; x A = \, x B = x c = \\ x A = ~, x A = \, x c = \\ and 
x A = xb — \i %c — 2> are developed to investigate the mixing enthalpies for ternary fee solid 
solution phases in the present work. The organization of this paper is as follows: the impact 
of ternary interaction parameters on a ternary solution phase in the CALPHAD approach 
is briefly reviewed. Then the generated ternary fee SQSs are characterized in terms of their 
atomic arrangements to reproduce the pair and multi-site correlation functions of completely 
random fee solid solutions. Finally, the generated SQSs are applied to the Ca-Sr-Yb system 
which presumably has fee solid solution phases throughout the entire composition range in 
all three binaries and ternary. 



II. TERNARY INTERACTION PARAMETERS 

The Gibbs energy of a ternary solution phase, (f>, is expressed as 6 

c c 

g 4, = +rtJ2 Xi ln Xi +xs ° hin ^ +xs Gtcrn,<t> w 

i=l i=l 

where X{ is the mole fraction of element i, is the Gibbs energy of phase of pure 
element i, xs Q hin ^ and xs G tcrn,4> are the excess Gibbs energies of the constitutive binary and 
ternary systems, respectively. The excess Gibbs energies for binary and ternary systems can 
be further described as 

c— 1 c n 

xs G Un = J2H •'•<•'•> E " ,: <M' - x tf ( 2 ) 

i=l j>i v=0 
c-2 c-1 c 

X s G tern = ^ ^ ^ X lXj X k { L U + L ^ + ^ ^ 
i=l j>i k>j 

where v Lf- is the v-th. order interaction parameter 13 in a binary and normally described as 

v Lt = v a + v bT (4) 
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where v a and v b are model parameters evaluated from experimental information. Ternary 
parameter of element i, Lf, in Eqn. 3 also have the form of Eqn. 4. If all three L-parameters 
in Eqn. 3 are identical, as in a ternary regular solution 6 , 

Lf = 4 = L* = L% k (5) 
then the ternary excess Gibbs energy shown in Eqn. 3 can be further simplified to 

^ tcrn = EEEw4 ( 6 ) 

i=l j>i k>j 

since Xi + Xj + Xk — 1 in a ternary. 

III. TERNARY FCC SPECIAL QUASIRANDOM STRUCTURES 

Thermodynamic properties of solid solution phases can be calculated in several ways 
(e.g., see Refs. 14 ' 15 ' 16,17 ' 18 ' 19 ). A popular approach is to use a database of first-principles cal- 
culations to determine the so-called Effective Cluster Interactions (ECI) that describe the 
energetics of the alloy system of interest. These interactions are then used as an input for ef- 
ficient statistical mechanics methods, such as the Cluster Variation Method or Monte Carlo 
simulations. While this general approach is able to model ordered phases (with potential 
point defects) and disordered (with potential short-range order) within a unified framework, 
it can be computationally demanding. Fortunately, in cases where the sole objective is 
obtain a reliable thermodynamic model for disordered phases that can be reasonably as- 
sumed to lack significant short-range order, the use of SQS provides a considerably more 
straightforward and computationally efficient approach. 

As discussed in Section I, first-principles study of SQS can effectively consider the limi- 
tations discussed above, and it has been demonstrated that binary SQSs can be applied to 
calculate thermodynamic properties of binary solid solutions, such as mixing enthalpy, for 
fee, bec, and hep structures 10 ' 11 ' 12 . It is thus anticipated that first-principles calculations of 
ternary SQSs should be able to reproduce thermodynamic properties of ternary solid solu- 
tions if their atomic configurations, which are represented as correlation functions, are very 
close to those of ternary solid solutions. Correlation functions of solid solution phases are 
well derived in Inden and Pitsch 20 . In the following section, the correlation functions for 
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binary and ternary systems are briefly summarized. 



A. Correlation functions 



The normalized correlation functions, Hfc, in crystalline structures are denned as 




2> 




■ ■ - a. 



.Cfc-l 

k 



(7) 



k site 



where the sum is over all the distinctive fc-site clusters, which are geometrically equivalent, 
in the N lattice sites structure. When fc=l, 2, 3, . . . , then fc-site clusters are point, pair, 
triplets, and so forth. Site occupation variables are denoted as a^, where the subscript fc 
indicates that the fc-th constituent located in the corresponding site. The superscript Ck takes 
values 2,3, • • • , C, with C as the number of constituents, which represents a constituent Ck 
at a given lattice site. 

For binary systems when C = 2, conventional values of the site occupation variables 
are ±1 depending on whether a lattice site is occupied by A or B atoms. According to Eqn. 
7, the normalized point correlation function for the second constituent site (B atom sites) 
is given as ITf = j^^cri^ with a± — 1 or -1. It is worth noting that the atom sites do 
not need to be distinguished in a binary system since they are switchable. With the mole 
fractions of A and B being xa and xb, respectively, (xa + xb = 1) then ITf = xa — %b- For a 
fc-site cluster, the normalized correlation functions for the binary solid solution is formulated 
as 



Similarly, for ternary systems when C — 3, the values of the site occupation variables Ok 
conventionally take +1, 0, or -1 if a lattice site is occupied by A, B, or C atoms, respectively. 
The normalized point correlation function for the second constituent site (B atom sites) is 
given as Tlf = -^J2 a i~ 1 with o\ = +1, 0, or, -1. For the third constituent site (C atom 
sites), the correlation function can be given as ilf = jj- of" 1 with <j\ = +1, 0, or, -1. With 
the mole fractions of A, B, and C being xa, %b, and xq, respectively (xa+xb+xc = 1); then 
Hi = xa~ %c and H\ = xa + xc- The vanishing of xb is due to its site occupation variable 
being 0. For a fc-site cluster with n# B atom sites and uq C atom sites (% + ric = fc), the 




k 



(8) 
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normalized correlation functions for the ternary solid solution is denoted as 

TT fc = (x A - x c ) nB {x A + x c ) nc with n B + n c = k (9) 
B. Generation of ternary fee SQS 

In the present work, two different ternary fee SQSs are generated. The first SQS is at the 
equimolar composition where xa = %b = %c = \ and the second \s at x a = \, x b = xc = \- 
By switching the occupation of the A atoms in the second SQS with either B or C atoms, 
two other SQSs can be obtained where xa = \, xb = \, xc = \ and x A — Xb — \, 
xq = \- Therefore, mixing enthalpy at four different compositions in a ternary system can 
be determined from first-principles total energy calculations of ternary fee SQSs by 

AH(A a B b C c ) « E{A a B h C c ) - x A E(A) - x B E(B) - x c E(C) (10) 

where E represents the total energy of each structure, and the reference states for all pure 
elements are fee. 

When the number of atoms in the ternary SQS is less than 24, the Alloy Theoretic 
Automated Toolkit (ATAT) 21 has been used to generate ternary fee SQSs. Since the ATAT 
enumerates all the atomic configurations within each supercell and then checks its correlation 
functions, the time needed to find SQSs increases exponentially as the size of a supercell 
increases. For the sake of efficiency, to find SQSs bigger than 24-lattice sites, a Monte-Carlo- 
like scheme 22 has been used. In each supercell with different lattice vectors, atom positions 
are randomly exchanged between the atoms and correlation functions of a supercell are 
calculated after every alternation. If the correlation functions of a new state are getting 
closer to those of random solutions, then the new configuration is accepted. Otherwise the 
new state is discarded and another configuration will be generated from the previous one. 
This process continues until the atomic arrangement of a supercell converges to the closest 
correlation functions of the random solution. In both methods, direct search via ATAT and 
Monte-Carlo-like scheme, a supercell whose correlation functions matches best with that of 
a completely random structure is chosen as the SQS at a given number of lattice sites. 

The selected SQSs at two different compositions, SQS-24 when xa = xb = xc = \ and 
SQS-32 when xa — \ and Xb — Xc — \, are shown in FIG. 1. These two SQSs are selected 
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for later calculations because they are adequate with respect to their size and correlation 
functions at each concentration 31 . The space groups of both structures are PI with all the 
atoms at their ideal fee sites. The correlation functions of the generated two SQSs are given 
in TABLEs I and II, respectively. 




(b) SQS-32 when xa — \, xb — %c = \ 

FIG. 1: Atomic arrangements of the ternary fee SQS structures in their ideal, unrelaxed forms. 
All the atoms are at their ideal fee sites, even though both structures have the space group PI. 

IV. FIRST-PRINCIPLES METHODOLOGY 

The Vienna Ab initio Simulation Package (VASP) 23 was used to perform the density 
functional theory (DFT) electronic structure calculations. The projector augmented wave 
(PAW) method 24 was chosen and the generalized gradient approximation (GGA) 25 was used 
to take into account exchange and correlation contributions to the Hamiltonian of the ion- 
electron system. An energy cutoff of 364 eV was used to calculate the electronic structures 
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TABLE I: Pair and multi-site correlation functions of ternary fee SQS structures when xa = ^b = 
xc = \- The number in the square bracket next to Hk, m is the number of equivalent at the same 
distance in the structure, the so-called degeneracy factor. 
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of all the SQSs. 5,000 Appoints per reciprocal atom based on the Monkhorst-Pack scheme 
for the Brillouin-zone sampling was used. 

In all first-principles calculations of binary and ternary SQSs in the present work, struc- 
tures are relaxed in two ways: full relaxation and volume relaxation to preserve the fee 
symmetry. It should be emphasized here that the preserved symmetry is that of SQS when 
all the atoms are substituted into a single element, which is the underlying fee symmetry. 
For the full relaxation of SQS, structures are completely relaxed with respect to all degrees 
of freedom i.e. cell shape, volume, and ionic positions, while the symmetry-preserving re- 
laxation only allows to change the cell volume for cubic structures, such as fee and bee. In 
calculating hep SQSs, however, relaxing only the volume will fix the c/a ratio of the under- 
lying hep symmetry 12 . Therefore, the shape of SQSs has to be relaxed as well as the volume 
in hep SQSs calculations. Since the present work focuses on the calculations of fee SQSs, 
symmetry-preserving relaxation is equivalent to volume relaxation. For symmetry-preserved 
calculations, all the atoms are still at the lattice sites of fee's and only the effective lattice 
parameter of fee changes with this constrained relaxation scheme. However, the local relax- 
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TABLE II: Pair and multi-site correlation functions of ternary fee SQS structures when xa = | ; 
%B = %c = \- The number in the square bracket next to Hk,m ls the number of equivalent at the 
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-0.01172 
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ation due to the like- and dislike-bondings cannot be taken into account. Further discussion 
with regard to the different relaxation schemes of SQSs can be found in Ref. 12 . 



V. RESULTS AND DISCUSSIONS 

The mixing enthalpy derived from first-principles calculations of ternary fee SQSs should 
be compared to relevant experimental measurements. However, there are always ordering 
effects at low temperatures where the fee solid solutions show complete solubility in binaries, 
such as Cu-Au and Au-Pd systems. Due to such ordering in binaries, it is almost impossible 
to reliably determine the mixing enthalpy for ternary fee solid solutions from experiments. 

In this work, the Ca-Sr-Yb system has been selected to apply the generated ternary fee 
SQSs which presumably has complete solubility in the fee phase for all binaries and ternary 
without any reported order /disorder transition. Both the Ca-Sr and Ca-Yb systems show 
complete solubility for both fee and bee phases at low and high temperatures respectively 
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without the formation of any intermetallic compounds 26,2 . There is no reported phase 
diagram for the Sr-Yb system, however, from the similarity of the two binary systems, 
Ca-Sr and Ca-Yb, it can be postulated that Sr-Yb also would have complete solubility for 
both fee and bee phases. Consequently, it may cautiously be expected that the combined 
ternary, the Ca-Sr- Yb system, would have the fee solid solution phase throughout the entire 
composition range at low temperatures. 

A. Binary SQSs for the Ca-Sr- Yb system 

Prior to applying the ternary SQSs to the Ca-Sr- Yb system, the mixing behavior of the 
fee phase in its constitutive three binaries was investigated through 8-atom binary fee SQSs 
at three different compositions, namely x=0.25, 0.5, and 0.75 in A\_ X B X alloys. Calculated 
mixing enthalpies from binary fee SQSs are combined with experimental data from the 
literature 28 ' 29 to evaluate interaction parameters in Eqn. 4 for each binary. For the sake of 
simplicity, parameters for bec have been modeled as identical to those of fee. Binary bec 
SQSs also have been calculated and included in the parameter evaluation process. FIG. 2 
shows that the formation energies of binary bec and fee SQSs very close to each other in all 
cases. The congruent melting of bec is observed in both the Ca-Sr and Yb-Ca systems, thus 
the Sr-Yb system has been modeled to have it as well on the assumption that Sr-Yb would 
behave similarly. 

Calculated mixing enthalpies for fee and bec phases of the three binaries are shown in 
FIG. 2 with first-principles calculations of binary fee and bec SQSs. Mixing enthalpy of the 
liquid phase in Ca-Sr is also calculated and compared with experimental measurement 29 . 
Structural analysis of binary SQS calculations shows that the local relaxation effect in both 
fee and bec SQSs is small after the full relaxation and the difference between fully relaxed 
and symmetry preserved structures are at most ~1 kJ/mol. 

It is intriguing to see that only the Yb-Ca system has the mixing enthalpy close to zero 
among three binaries in FIG. 2, which implies that Yb-Ca is likely to have ideal mixing in 
both fee and bee phases. The ideal mixing of Yb and Ca is attributed to their similar lattice 
parameters. The other two systems, Ca-Sr and Sr-Yb, have rather positive (~2 kJ/mol) 
mixing enthalpies and significant lattice parameter differences (ac a ~ «Yb <C ctsr); which can 
be the reason for the non-ideal mixing behavior, unlike Yb-Ca. 
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(c) Sr-Yb (fee) 



(d) Yb-Ca (fee) 



FIG. 2: Calculated mixing enthalpies for (a) the liquid phase in the Ca-Sr with experimental 
measurement 29 and the fee phase for (b) Ca-Sr, (c) Sr-Yb, and (d) Yb-Ca with first-principles 
calculations of fully relaxed SQSs. Note that fee and bec phases are evaluated with same model 
parameters. Open and closed symbols represent symmetry preserved and fully relaxed calculations, 
respectively. 

The bond length analysis for the fully relaxed SQSs in FIG. 3 shows that first nearest- 
neighbor average bond lengths follow Vegard's law closely in all calculations. This obser- 
vation indicates that the lattice parameter of the fee solid solution varies linearly with the 
composition change and there is no significant geometrical distortion. 
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FIG. 3: SQS bond lengths for three binaries, Ca-Sr, Yb-Sr, and Yb-Ca. Error bars correspond to 
the standard deviation of the bond length distributions. 

B. Ternary fee SQSs for the Ca-Sr- Yb system 



First-principles calculations of ternary fee SQSs at four different compositions in the 
Ca-Sr- Yb system, namely xq 



3 ) -^Ca 



^Sr 



Xyb 



4> x Sr 
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2 • 



xca, = xyb = \\ and xy^ = \i Xca, = ^Sr = \i have been considered to investigate the 
ternary interactions. Three isoplethal sections, connecting the equimolar composition and 
three other compositions when = 1/2, Xj/x^ = 1, are selected to examine the mixing 
enthalpy for the Ca-Sr- Yb ternary system. 

Calculated enthalpies of mixing from ternary fee SQSs are shown in FIG. 4, including 
extrapolated results from the three binaries and improved enthalpies of mixing to reproduce 
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FIG. 4: Calculated mixing enthalpies for the fee phase in the Ca-Sr-Yb system with first-principles 
calculations results of ternary SQSs in three pseudo-binaries, Ca-SrYb, Sr-YbCa, and Yb-CaSr. 
Extrapolated mixing enthalpy for the fee phase from three binaries is compared with the ones with 
ternary interaction parameters. Open and closed symbols represent symmetry preserved and fully 
relaxed SQS calculations, respectively. 



ternary fee SQSs results by introducing ternary interaction parameters. All the fully relaxed 
ternary SQSs show that the effect of local relaxation is also small as in its constitutive 
binaries, thus the energy differences between the symmetry preserved and fully relaxed 
calculations are quite small. FIG. 5 shows the radial distribution analysis of the fully relaxed 
SQS at the equimolar composition. The narrow distribution along each of the bond-lengths 
indicates that the effect of local relaxation is small. As shown in FIG. 4, energy differences 
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between fully relaxed and symmetry preserved calculations are small(~ 0.5 kJ/mol). 
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FIG. 5: Radial distribution analysis of CaiSriYbi ternary fee SQSs. The dotted line in (b) under 
the smoothed and fitted curves are the error between the two curves. 

As can be seen in FIG. 4, the mixing enthalpy extrapolated from the three binaries 
is slightly lower than that derived from first-principles calculations of ternary fee SQSs. 
Thus, ternary interaction parameters are introduced to improve the ternary mixing enthalpy. 
According to Eqn. 3, the contribution from the ternary excess Gibbs energy for the fee phase 
in the Ca-Sr-Yb system can be denoted as 
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X s G temM = XCaXSrXYb (L^X Ca + L^X Sl + Lf b X Yb ) (11) 

or 



xs srtern,fcc _ rfee (TO\ 

<J — ^Ca^Sr^Yb-^CaSrYb 

as simplified in Eqn. 6 when the ternary fee is considered as a regular solution. When three 
independent ternary interaction parameters (Lea = 25940, Ls r = 2913, and Lyh = —8645 
J/mol) are used, slightly better agreement with ternary SQSs was made than with a single 
regular interaction parameter (I/CaSrYb = 6736 J/mol). The calculated mixing enthalpy at 
the equimolar composition are evaluated as the same value regardless of the interaction 
parameters since all the data are equally weighted. 
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FIG. 6: First nearest-neighbor average bond lengths for the fully relaxed ternary fee SQSs in the 
Ca-Sr-Yb system. Error bar corresponds to the standard deviation of the bond length distributions. 

Six different bond lengths, three like-bondings and three dislike-bondings, of four fully 
relaxed SQSs have been analyzed. In FIG. 6 the bond lengths corresponding to the first 
nearest neighbors for all SQSs are presented. In all SQS calculations, the Sr-Sr bonding is 
always the longest and this is attributed to the biggest lattice parameter of Sr among the 
three elements. 

As a final analysis of the ability of the generated SQSs to reproduce the properties of 
ternary fee solid solutions, FIG. 7 shows the alloying effects on the electronic DOS in the Ca- 
Sr-Yb system. Three pseudo-binaries, connecting a pure element and equimolar composition 
of the other two elements, have been selected. Since all three elements have the same number 
of valence electrons, significant changes in the electronic DOS were not observed. Instead, 
the electronic DOS of the outer states gradually transformed into that of the pure element 
as the composition changed towards the pure element. The peak around ~ 9 eV of SrYb 
in FIG. 7(a) flattens out as Ca increases, and the peaks around ~ 8 eV and ~ 7 eV are 
become more pronounced as the content of Sr and Yb increases as shown in FIGs. 7(b) and 
7(c), respectively. 
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FIG. 7: Calculated density of states for Ca-Sr-Yb as three pseudo-binaries. 

VI. CONCLUSION 

In the present work, two ternary fee SQSs at different compositions, xa = %b = %c — \ 
and xa = \, xb = xc = \, are generated and their correlation functions are satisfactorily 
close to those of random fee solid solutions. The generated SQSs are applied to the Ca-Sr-Yb 
system which presumably has a complete solubility range without order/disorder transitions 
in ternary fee solid solutions. Mixing enthalpies for the fee phase in three binaries are 
evaluated from first-principles calculations of fee and bee SQSs and available experimental 
data with the CALPHAD approach. It is found that the local relaxation effect of fee and 
bec phases are very small and mixing enthalpies are slightly positive in all cases. Evaluated 
mixing enthalpies for the fee phase in three binaries are then extrapolated to the ternary 
system. 

First-principles results of four ternary SQSs at xc& = x$ r = xyh = §; ^Ca = \, ^Sr = 
x Y b = \\ £sr = f, £ C a = a^Yb = \\ and x Yh = §, «Ca = ^Sr = \ show that the local 
relaxation effect is also very small in the ternary system. Extrapolated mixing enthalpy 
from its constitutive binaries are slightly lower than those from first-principles calculations 
of ternary fee SQSs. Thus ternary interaction parameters for fee solid solution phases are 
introduced to further improve the ternary mixing enthalpy. It can be concluded that the 
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generated ternary fee SQSs are able to reproduce thermodynamic properties of ternary fee 
solid solutions and can readily be applied to other systems. 
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APPENDIX A: EVALUATED INTERACTION PARAMETERS 

Interaction parameters for the liquid and solid solution phases (fee and bee) in the three 
binaries of the Ca-Sr-Yb system — as evaluated from the CALPHAD modeling in the present 
work — are listed in TABLE III. Notations are explained in Eqn. 2. 



TABLE III: Evaluated binary interaction parameters for the Ca-Sr-Yb system (all in S.I. units). 
Gibbs energies for pure elements are from the SGTE pure element database 30 . 



Systems 


Phases 


Evaluated parameters 


Ca-Sr 


Liquid 


°L = 1680 + 2.68T, 1 L = 388 - 1.153T, 2 L = -856 + .631T 




fcc,bcc 


°L = 6511, 1 L = -382 


Sr-Yb 


Liquid 


°L = 8850 




fcc,bcc 


°L = 9624, 1 L = -488 


Yb-Ca 


Liquid 


°L = 2689, 1 L = 676 




fcc,bcc 


°L = 3207, 1 L = 280 



APPENDIX B: SPECIAL QUASIRANDOM STRUCTURES FOR THE 
TERNARY FCC SOLUTION PHASE 

Special quasirandom structures are Y-atom per cell periodic structures designed to have 
correlation functions close to those of completely random alloys. The ternary fee SQSs used 
in this work are presented. Lattice vectors are given as a, b, and c, and atom positions of 
A, B, and C atoms are given as A i: B^ and Cj, respectively. 



17 



TABLE IV: Structural descriptions of the SQSs for the ternary fee solid solution. Lattice vectors 
and atomic positions are given in fractional coordinates of the fee lattice. Atomic positions are 
given for the ideal, unrelaxed fee sites. 



A1B1C1 (SQS-24) A2B1C1 (SQS-32) 
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